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If supersymmetry is discovered at the LHC, the next question will be the determination of the 
underlying model. While this may be challenging or even intractable, a more optimistic question is 
whether we can understand the main contours of any particular paradigm of the mediation of su- 
persymmetry breaking. The determination of superpartner masses through endpoint measurements 
of kinematic observables arising from cascade decays is a powerful diagnostic tool. In particular, 
the determination of the gaugino sector has the potential to discriminate between certain mediation 
schemes (not all schemes, and not between different UV realizations of a given scheme). We recon- 
struct gaugino masses, choosing a model where anomaly contributions to supersymmetry breaking 
are important (KKLT compactification), and find the gaugino unification scale. Moreover, recon- 
struction of other superpartner masses allows us to solve for the parameters defining the UV model. 
The analysis is performed in the stop and stau coannihilation regions where the lightest neutralinos 
are mainly gauginos, to additionally satisfy dark matter constraints. We thus develop observables to 
determine stau and stop masses to verify that the coannihilation mechanism is indeed operational, 
and solve for the relic density. 



I. INTRODUCTION 

The Large Hadron Collider (LHC) will soon be test- 
ing many ideas for physics beyond the Standard Model 
(SM). Of these, low energy supersymmetry [T] is the best 
motivated candidate for new TeV scale physics. 

If supersymmetric partners are indeed detected at the 
LHC, the next step would be understanding the under- 
lying scheme by which supersymmetry breaking is medi- 
ated to the visible sector. Of course, the full determina- 
tion of an exact model of supersymmetry breaking and 
mediation may be challenging, but one can hope that the 
broad contours of the scheme can be understood. 

How soon (if at all) will we be able to distinguish the 
underlying mechanism by which supersymmetry break- 
ing is mediated to the visible sector? 

It is useful to first understand the above question from 
the perspective of distinguishability studies undertaken 
from a top-down point of view |2] • Broadly, the program 
is to fit a set of UV parameters defining a sufficiently 
tractable model of supersymmetry breaking and media- 
tion with observables measured in the experiment, and 
ask questions about distinguishability in the space of UV 
model parameters. Some models one can consider in such 
distinguishability tests are (i) mSUGRA or CMSSM [3] 
(defined at the GUT scale by universal scalar and gaug- 
ino masses and trilinear couplings mo,mi/2, and Aq, as 
well as tan/3 and ratio of Higgs vevs), (iz) the minimal 
anomaly mediated supersymmetry model [1] (defined by 
mo, tan/3, and universal superpartner mass scale maux), 
(iii) minimal gauge mediation (defined by the messenger 
scale Mmess, supersymmetry breaking scale A and tan/3), 
as well as some examples inspired by string phenomenol- 



ogy: (iv) Large Volume Scenarios j5j (defined by tuq and 
tan/3), and (v) mirage mediation models [B] in the context 
of KKLT [7]. 

However, the important point is that a given set of 
measured observables will in general not only point back 
to different UV models, it may point to different media- 
tion schemes. It may be more optimistic to ask whether 
we can at least distinguish between more gross proper- 
ties, such as mediation schemes. 

It is important, then, to first discern whether a partic- 
ular scheme of supersymmetry breaking mediation (such 
as, say, anomaly contributions) shows up unmistakably 
in the low energy spectrum, before locking into a par- 
ticular UV model and undertaking a distinguishability 
study in the UV model parameter space or attempting 
to constrain it. 

Which masses should one attempt to solve to this end? 

One can look for relatively model-independent soft 
mass patterns that hold out clues for mediation schemes. 
As emphasised in |8] , the gaugino sector offers a particu- 
larly promising portal for understanding the underlying 
mechanism of mediation. This is mainly because within 
the context of the Minimal Supersymmetric extension of 
the Standard Model (MSSM), the quantity Ma/gl does 
not run at one loop (where a = 3, 2, 1 refers to the SM 
gauge groups SU{3) x SU{2) x C/(l), while Ma refer 
to gaugino masses and ga to gauge couplings). If one 
makes two further assumptions: (i) that gaugino masses 
are dominated by contributions determined by tree level 
gauge kinetic functions and (ii) these tree level contri- 
butions are universal, then one obtains the well known 
mSUGRA pattern, with gaugino unification at the GUT 
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scale 



mSUGRA : Mi : Nh : M3 



1:2:6 



(1) 



This pattern is obtained in a variety of schemes: (i) grav- 
ity mediation, with different UV scenarios (m) gauge me- 
diation (iii) gaugino mediation [9]. We wiU elaborate 
further in the next Section. 

In a previous study [T^, we have described a series 
of measurements of superpartner masses at the LHC for 
mSUGRA models using end-point techniques [Tl . 

On the other hand, gaugino masses may be domi- 
nated entirely by one loop contributions coming from 
the SUGRA compensator, determined by the confor- 
mal anomaly of the effective theory at TeV scale. At 
some level, the flavor problem embedded in dominant 
tree level contributions is an indication that tree level is 
perhaps not the end of the story. For pure anomaly me- 
diation, Ma{Q) / g'^iQ) is proportional to the respective 
beta- function coefhcients at the scale Q. Anomaly me- 
diation effects are always present and become dominant 
over gravity mediation when the supersymmetry break- 
ing sector is sequestered from the visible sector; however, 
the scheme by itself suffers from tachyonic sleptons. 

Mirage mediation is a hybrid of the mSUGRA pattern 
and anomaly pattern, in which the tree level contribution 
and the one loop contribution are both equally compet- 
itive. The quantity Ma/g^ then depends both on the 
universal mSUGRA contribution, as well as the RG beta 
coefficients, with comparative strengths given by a pa- 
rameter a: 
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with AIq ~ 1 TeV a mass scale. 

The above facts translate into the following ratio of 
low energy gaugino masses 

Anomaly : Mi : M2 : M3 - 3.3 : 1 : 9 , 
Mirage : Mi : M2 : M3 

(1 + 0.66a) : (2 + 0.2a) : (6 - 1.8a) (3) 

Clearly, a deviation of the measured gaugino spectrum 
from the mSUGRA pattern would signal the relative im- 
portance of one loop contributions (which, of course, 
could have a plethora of model origins). 

A. Mass Measurements at the LHC 

In this paper, we take the first bottom-up steps toward 
the diagnosis of mediation schemes at the LHC, by choos- 
ing a representative model and attempting to reconstruct 
low energy masses. 

Mass reconstruction gives us the following information: 
(z) As we have outlined above and spell out in more 
detail later, mass measurements can distinguish between 
mediation schemes. Of course, the choice of the UV 



model must be judicious - for example, reconstructing 
just the gaugino sector would still leave one incapable of 
distinguishing between the various models and schemes 
which give rise to the mSUGRA pattern of masses (this 
has also been verified by distinguishability tests [H]). We 
thus choose a model (KKLT compactification) where the 
imprints of anomaly mediation are unmistakable. 

One important point in this paper is that in cer- 
tain regions of parameter space, the crucial information 
from the gaugino spectrum regarding the scheme may 
be gleaned at relatively lower luminosity by a judicious 
choice of observables. Meanwhile, the details of a par- 
ticular model only become apparent as other masses are 
reconstructed. 

(a) The measurement of low energy masses enables 
us to make statements about the dark matter relic den- 
sity. In models with i?— parity invariance, the lightest 
superpartner (typically the lightest neutralino xi°) is a 
dark matter candidate, and one would like to compute 
the relic density ilh^ associated with it. The direct mea- 
surement of stop (t) and stau (f) masses enables us to 
verify whether we are in regions of parameter space where 
stop-neutralino {i — xi°) and stau-neutralino (f — xi'^) 
coannihilation effects are important. We note that the in- 
vestigation of mediation schemes, the first goal described 
above, typically allows one to be in much larger areas of 
parameter space, and one need not necessarily work in 
coannihilation regions. However, to satisfy relic density 
constraints, they are crucial. 

Previously, we have undertaken studies of mass mea- 
surements in the t — xi° coannihilation region in the 
context of the mSUGRA model [lOj il3j . The measure- 
ment of third generation squark masses presents its own 
challenges. Reconstruction of stop and sbotttom (5) is 
very hard in a cascade decay chain since both stops and 
sbottoms decay into b quarks. Further, to make the sit- 
uation worse, in the stop coannihilation region the stop 
decay produces a lower pT jet due to the proximity of 
the lighter stop and the lightest neutralino masses. We 
invoke two new observables to measure lighter stop and 
sbottom masses in the cascade decays. 

Thus, in this paper we fulfil the following objectives: 

(i) Constructing observables using different combina- 
tions of the jets, r's and W's in the final states to 
determine the masses of supersymmetric particles in a 
model where anomaly contributions are important. We 
construct observables in the stop-neutralino and stau- 
neutralino coannihilation regions. 

(a) Using the observables, we solve for the masses of 
the gluino, the two lighter neutralinos, the squarks (of 
the first two generation) and the lightest stau. We also 
determine the lighter stop and sbottom masses in the 
stop coannihilation region. 

(iii) Using the neutralino and gluino masses we deter- 
mine the gaugino unification scale, thereby establishing 
the imprint of one-loop anomaly contributions competi- 
tive with tree level contributions (mirage pattern) from 
experimentally measurable observables. 
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(iv) Using the masses of the other particles in conjunc- 
tion with the gauginos, we (a) determine the parameters 
defining a particular UV completion of the mirage pat- 
tern (KKLT compactification with visible sector on D7 
branes) and (b) test whether we are in a coanihilation 
region and dark matter relic density is satisfied. 

A caveat is in order. The observables we construct in 
the current work enable us to measure the masses of the 
gluino and the two lighter neutralinos, which should be 
mainly gauginos for us to make statements on mediation 
within the structure outlined above. In the coannihila- 
tion regions, this is guaranteed, and thus we are able to 
simultaneously solve the relic density as well as deter- 
mine the mediation scheme. It would be very interest- 
ing to probe similar issues in regions of parameter space 
where the Higgsino component in the neutralinos is not 
insignificant. However, such a study is beyond the scope 
of this paper. 

Another comment pertains to the fact that we will be 
solving the model parameters from the masses. As we 
have stressed, our aim in this paper is to show that the 
determination of a given set of masses can unmistakably 
establish the contours of a mediation scheme (point {in) 
above), while it may point back to different models. In 
point {iv) above, we will thus not attempt a distinguisha- 
bility test of our particular model (KKLT) - rather, we 
will fit our model parameters with the masses we solve 
from the observables. 

The plan of the paper is as follows. In Section |llj we 
describe the contributions to gaugino masses, describe 
our choice of the UV model in more detail, and fix our 
benchmark points. In Section III we develop the kine- 
matic observables to measure sparticle masses in the stop 
coannihilation region. Included in this section is the de- 
velopment of specific observables required for the recon- 
struction of stop and sbottom masses. In Section IV we 



give results for the gaugino mass pattern and unification 
scale and solve for the other masses and the full parame- 
ter set defining the UV model, along with the relic den- 
sity. Moreover, we solve for the stop and sbottom masses. 
In Section |Vj we perform the above analysis in the stau 
coannihilation region of parameter space. We end with 
our conclusions. 



A. Contributions to Gaugino Masses 

The various contributions to Maj g\ 4D effective su- 
pergravity are 

— = Ma ^Ma lanomaly+Ma | other, (4) 

where m}'^'^ = 1/2F^ diff'\ with signifying the non- 
zero F— component of a hidden sector field Xj breaking 
supersymmetry, and fa denoting the visible sector gauge 
kinetic functions. The one loop anomaly term receives 
contributions from the conformal and Konishi anoma- 
lies, with the conformal anomaly contribution given by 

conformal — 

(l/167r2)fe^(F^/C), with C denoting 
the chiral compensator. We will not deal with the Kon- 
ishi anomaly contribution in any detail in this paper. 

The other one loop contributions Ma^ ^ I other encapsulate 
field theoretic gauge threshold and UV-sensitive contri- 
butions like KK thresholds to the gaugino masses. We 
will assume that UV-sensitive contributions are subdom- 
inant (the assumption is purely because they are more 
model-dependent) . 

The mSUGRA pattern, Eq. [T] arises in situations 
where Malg\ dominated by universal tree level piece 

Ma ■ This happens in a variety of schemes: (i) gaugino 
mediation in higher dimensional brane models (ii) grav- 
ity mediation, with different UV scenarios such as dila- 
ton/moduli mediation in heterotic string theory or large 
volume compactification in type IIB tring theory. Gauge 
mediation also gives the mSUGRA pattern of gaugino 
masses, since the gaugino masses are dominated by uni- 
versal one-loop gauge threshold contributions. 

On the other hand, the mirage pattern in which tree 
level and loop level contributions are competitive occurs 
in gravity mediation in various UV scenarios, for exam- 
ple KKLT flux compactification with visible sector on 
Dl branes and explicit supersymmetry breaking by anti- 
i'S branes or spontaneous breaking by a matter sector. 
Mirage patterns are also obtained in deflected anomaly 
mediation [H] . 

However, the important point is that deviations from 
the mSUGRA pattern may signal some anomaly contri- 
bution at work. 

We will take as an example the UV completion of 
KKLT with visible sector on Dl branes and explicit su- 
persymmetry breaking by an anti-_D3 brane. 



II. MODEL AND BENCHMARK POINTS 

As emphasised in Section [ij a given pattern of low 
energy gaugino masses accommodates diverse mediation 
schemes, not to mention multiple models. We mention 
below a catalogue of such models and schemes, outlined 
in [H], to better situate our particular UV completion. 



B. KKLT with D1 branes: The UV Model 
Parameters 

In this subsection, we sketch the outlines of the UV 
model, only with a view to defining the model parameters 
we will be working with. For details, we refer to [7j. 

The basic elements in a KKLT-type model of string 
compactification are: (i) background fluxes on a type 
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IIB Calabi-Yau three-fold giving a Gukov-Vafa-Witten 
superpotential contribution, and (ii) gaugino condensa- 
tion on D7-branes or Euclidean _D3 instantons giving a 
non-perturbative superpotential contribution. These sta- 
bilize complex structure moduli and the dilaton, as well 
as Kahler moduli, in an AdS vacuum. Supersymmetry 
breaking by an anti-I?3-brane then lifts the solution to a 
de Sitter vacuum. The visible sector is constructed with 
D7 branes in the bulk. 

The soft masses at the GUT scale are 



Ma = Mo+'^bagl 



167r2 

"•^3/2 
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where Mq and rhi are pure tree level modulus contribu- 
tions, given as functions of the Kahler modulus T, while 
other terms are one- loop anomaly contributions. 1713/2 is 
the gravitino mass. In the above, 



ba = -3tr(r2(Adj))+5]tr(r,2(</.,)) 

i 

a jk 

n 2 dj, 

a jk 



(6) 



where the quadratic Casimir Cf((/)i) — {N'^ — 1)/2A^ 
for a fundamental representation (pi of the gauge group 
SU{N), C!^{(j){) = qf for the U{1) charge q, of and 
J2ki VikiVjki is assumed to be diagonal. 

To list the set of parameters that give the low energy 
mass spectrum, we define the ratios 



C. Benchmark Points 

The parameter space of the KKLT model has been 
studied in [15], and we choose points that satisfy the 
stop-neutralino and stau-neutralino coannihilation con- 
straints, which appear in ample parts of model parameter 
space. We use darkSUSY jl6J to select exact benchmark 
points in the above regions. 

The benchmark point for the stop coannihilation re- 
gion is shown in Tables [ij [llj 



TABLE L Model parameters chosen at a stop coannihilation 
benchmark point. All masses are in GeV. 



Parameter Value 



a 
"^3/2 

tan/3 



4.5 
14000 
0.0 
0.5 
30 



TABLE IL Spectrum at a stop coannihilation benchmark 
point. All masses are in GeV. 



Particle 


Mass 


Particle 


Mass 


Particle 


Mass 


dL 


653.13 


Sl 


436.75 


X? 


286.21 


da 


635.86 


eR 


411.28 


X§ 


338.21 


UL 


647.91 


Tl 


315.08 


xg 


477.35 


Ur 


634.96 




417.70 


xl 


502.68 


h 


520.46 






xt 


337.32 




596.25 






xt 


500.41 


ii 


338.55 






9 


649.78 




616.22 











"^3/2 



MoHMpi/my^)'' 



1 — 71,: — 



Ml 



(7) 



where a represents the anomaly to modulus mediation 
ratio, while rii are modular weights, that parameterize 
the pattern of the pure modulus mediated soft masses. 

For convenience, we will henceforth choose to rescale 
a as follows 
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c forth 



ln(Mp/m3/2) a 



(8) 



From Eq. [5] then, it is clear that the parameters defin- 
ing the low energy soft masses of the model are 



nT-3/2, a, "m, riH, tan/3. 



(9) 



where we choose 7713/2 in place of Afo, split the matter 
modular weights into a universal sfermion weight rim and 
a Higgs weight uh- We have also replaced the Higgs mass 
parameters /i and B by tan/3 and Mz- 



The benchmark point for the stau coannihilation re- 



gion is shown in Table III 



TABLE III: Model parameters chosen at a stau coannihilation 
benchmark point. All masses are in GeV. 



Parameter Value 



a 
"13/2 

riH 
tan/3 



7.5 
10000 
0.5 
1.0 
30 



The spectrum at the stau coannihilation benchmark 



point is shown in Table IV 



Note that our methods are valid in general, and the 
above benchmark points will be explored as an illustra- 
tion. In particular, we will also study benchmark points 
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TABLE IV: Spectrum at a stau coannihilation benchmark 
point. All masses are in GeV. 



Particle Mass Particle Mass Particle Mass 





845.49 


eL 


426.91 


X? 


284.17 


dR 


813.52 


e-R 


367.70 


xl 


389.17 


Ul 


841.39 


fi 


309.75 


xi 


548.88 


Ur 


815.27 




425.68 


xl 


569.04 


61 


735.87 






xt 


389.32 


62 


791.30 






xt 


568.10 


ii 


600.23 






9 


897.55 


i2 


810.20 











with higher gluino mass, preferred by current LHC data, 
in Sections |IV C| and |V C[ where we will show that we 
need larger luminosity to establish the same set of ob- 
servables. 

Since tan/3 is on the large side for the benchmark 
points, the lighter stau mass is between the lightest and 
next to Hghtest neutralinos. The lightest neutralino is 
mostly Bino and the next to lightest is mostly Wino. The 
Higgsino components are negligible since the dark matter 
relic density is satisfied by the coannihilation mechanism. 

III. KINEMATIC OBSERVABLES IN THE STOP 
COANNIHILATION REGION 

In this section, we present the measurement of physi- 
cal observables that will be used to solve for the gaugino 
masses, as well as the masses of the f and q. Moreover, 
we construct observables that will be needed to solve for 
the t and h masses. In the next section, we will present 
the mass and model parameter solutions, and also present 
gaugino unification. We will also give results for a bench- 
mark point with heavier gluino there. 

The mass spectrum of the model is determined using 
ISASUGRA [13. The spectrum is then fed to PYTHIA [T5] . 
which generates the Monte Carlo hard scattering events 
and hadron cascade. These events are passed to the de- 
tector simulator PGS4 HH]. 

A. 2 Jets + 2t + fh 

At the LHC, the main production processes for this 
model are gg, gq and qq. 

The relevant decay chains are: 

9 — > QlQl 

and 

QL X29 r^T^q XiT^T^q ■ (10) 

There is a high px jet, missing transverse energy, and a 
pair of oppositely charged r leptons. 



The following cuts were set to select events for this 
signal: 

(i) Missing transverse energy fJr > 180 GeV; 

la) PT.jetl + PT.jet2 |^ ^ 600 GcV; 

{in) Leading jet cuts: At least two jets should be 
present, each with pr > 200 GeV in \r]\ < 2.5. The 
jets are both required to be non ^-tagged jets, and the 
event is discarded when either of them is tagged as a 6 
jet; 

(iv) Soft jet cuts: Any jet with px > 30 GeV in I77I < 
2.5 is accepted in the analysis; 

(v) T cuts [5D]: At least two t leptons, with visible 
PT > 15 GeV in \r]\ < 2.5; 

Various kinematic observables can be constructed in 
this signal. Of them, we choose four that give the most 
precise endpoint measurements with 50 fb~^: M™^, 
M™'^, and the pr distributions of the higher and lower 
energy r's. 

The fifth observable (required to solve the five model 
parameters) is the peak of the Mcs distribution. 

1. M^'f 

The main challenge in measuring the invariant di-tau 
mass is the background created by uncorrelated r pairs. 
From the decay chain in question, it is clear that sim- 
ilarly charged r's are definitely uncorrelated and hence 
model this background quite well. Thus, we perform the 
opposite-sign (OS) minus like-sign (LS) subtraction on 
r's to obtain the signal. 

Figure [l] shows the histogram graph we obtained for 
the Mrr distribution. The luminosity is 50 fb~^, and we 
can see that a clear end-point is obtained. 




20 40 60 80 100 120 140 



Mx. (GeV) 

FIG. 1: Distribution of Mtt at a stop coannihilation bench- 
mark point. The pink (grey) histogram is composed of OS 
T pairs. The filled dot-dashed pink (grey) histogram is com- 
posed of LS r pairs and is normalised to the shape of the 
long tail of the the pink (grey) OS histogram. The OS— LS 
subtraction produces the black subtracted histogram. This 
subtracted histogram is then fitted with a straight line to find 
the endpoint of the distribution. The result for the endpoint 
is 49.71 ± 0.2(Stat.) GeV. The luminosity is 50 fb"^ 
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2. M;;° 

The T leptons from each event are sorted into OS and 
LS pairs. Then, each pair is combined with the leading 
jets from the same event, and the invariant mass distri- 
bution Af|™® is obtained. At this stage the Bi-Event 
Subtraction Technique (BEST) pT] is initiated, by first 
combining the r pair with leading jets from a separate 
event (which we also call a bi-event), to form the dis- 
tribution Mj'^""™"*. Since jets from a separate event 
are kinematically uncorrelated with the r pair from the 
current event, this bi-event distribution models the jet 
background very well. 

We note that the greater the number of different-event 
jets included in the bi-event distribution, the better the 
modelled background becomes. We include sufficient 
number of different-event jets to ensure a smooth back- 
ground distribution. 

After this, two subtractions are performed to obtain 
the final signal: (i) The OS— LS subtraction, which gets 
rid of uncorrelated r pairs and (ii) The BEST subtrac- 
tion, which gets rid of uncorrelated jets in the back- 
ground. In Figure [2j we show the Mj^r distribution 
where we find a very clear end point after the BEST 
subtraction. 




200 400 600 800 



Mj,, (GeV) 

FIG. 2: Distribution of M^tt at a stop coannihilation bench- 
mark point. The green (grey) same-event histogram is con- 
structed by combining each OS— LS t pair with a leading jet 
from the same event. The filled dot-dashed green (grey) bi- 
event histogram is obtained by combining the OS— LS r pair 
with a leading jet from a difi^erent event and is normalised to 
the long tail of the same-event histogram. The same-event 
minus bi-event subtraction (BEST) produces the black sub- 
tracted histogram. The subtracted histogram is fitted with a 
straight line to obtain the endpoint. The result for the end- 
point is 269.09 ± 3.18(Stat.) GeV. The luminosity is 50 fb'^ 



the pt of the t's in the ditau pair. The slope of the visi- 
ble pt distribution of the lower energy visible r from the 
di-tau pair (which we will denote by slope{pT)) is gen- 
erally a good observable, that carries information about 
the masses of f and Xi (in the case that the lower energy 
r is from f — > Xi+ t). If the stau mass is very close to 
the neutralino mass (in the case of stau-neutralino coan- 
nihilation) we have a low energy t and the slope of the 
Pt distribution is a good observable. 

Conversely, if the visible pt of a r is large, the slope 
of the Pt distribution is not a good observable. This is 
because the slope in the case of a higher energy r doesn't 
show enough variation. In such cases, we use the kine- 
matic information of the high energy r by combining it 
with a leading jet to construct Mjr- We will use the 
endpoint of the Mj^ distribution as an observable in the 
stau-neutralino coannihilation case, where one of the r's 
has large pt compared to the other. 

If the stau is somewhat in between the two neutralinos, 
the slopes of the pt distributions of both taus become 
good observables (provided, of course, that neither px is 
too large) as happens in this case. In such a scenario, we 
get two observables depending on two neutralinos and 
the stau mass without involving any jet. It is convenient 
to define the new observables Pt, am and PT.diffi which are 
the arithmetic mean and difference of the slopes of the 
higher and lower px t's. These variables can typically be 
measured down to lower luminosity than Mjr, since they 
do not involve the subtractions required for observables 
involving jets, as we describe below. 

To obtain the distributions, first the OS— LS subtrac- 
tion is performed. Then, we take the logarithm of the 
PT,high and pt,Iow distributions for the higher and lower 
energy r respectively. The histograms are fitted using 
a linear function to find the slopes. Finally, the slopes 
are combined into the arithmetic mean pt,am and the 
difference pT.difF^ 

Pt,am = ^ {slope{pTMgh) + slope{pT.\ow)) 

PT.difi = ^ (s'ope(pT,high) - slope{pT,iow)) ■ (H) 

Note that these observables, formed by combining the pt 
information of the two r's, are particularly necessary in 
the case when they are close in mass, since we do not 
know if a given r is originating from the f decay or the 
X2 decay. 

In Figure [3] and Figure |4j we show the pr distributions 
of two r's present in the sample. 



B. 4 jets + pT 

3. Slope ofpT ofT r^Yie Moff distribution is formed by combining the pr 

of the first four leading jets and the missing energy 

The mass of the f lies between the two lightest neu- 
tralinos, and its exact location approximately determines Mcs — PTjeti + PT,jet2 + PT,jet3 + PT,jet4 + fr (12) 
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FIG. 3: Distribution of pT.high at a stop coannihilation bench- 
mark point. The pink (grey) histogram is composed of OS r 
pairs. The filled dot-dashed pink (grey) histogram is com- 
posed of LS r pairs and is normalised to the shape of the 
long tail of the OS histogram. The OS— LS subtraction pro- 
duces the black subtracted histogram. This subtracted his- 
togram is then fitted with a straight line to find the slope of 
the distribution. The value of the slope{pT, high) observable is 
-0.0522 ± 0.0021(Stat.). The luminosity is 50 Ib'^ 




Pt,,ow 

FIG. 4: Distribution of Pt.Iow at a stop coannihilation bench- 
mark point. The pink (grey) histogram is composed of OS 
T pairs. The filled dot-dashed pink (grey) histogram is com- 
posed of LS r pairs and is normalised to the shape of the 
tail of the OS histogram. The OS— LS subtraction produces 
the black subtracted histogram. This subtracted histogram 
is then fitted with a straight line to find the slope of the 
distribution. The value of the slope{pTAo-w) observable is 
-0.1178 ± 0.0040(Stat.). The luminosity is 50 fb-\ 

These jets effectively come from the ghiino and squark 
decays. 

The effective mass Mcs is constructed from the 4jets 
-|- ^ sample, with the following cuts 

(i) Number of jets: At least four jets in the event; 

{a) Leading jet cuts: The first two leading jets each 
have pT > 200 GeV in \r]\ < 2.5; 

(Hi) Soft jet cuts: Jets with pT > 30 GeV in |?7| < 2.5 
are accepted in the analysis; 

{iv) The jets are not b tagged; 

(v) fk > 180 GeV; 



(vi) PT.jeti + PT.jet2 + fk > QOO GeV; 

(vii) No e's and /i's with px > 15 GeV; 
(via) Transverse sphericity St ^ 0.2 . 

In Figure [5j we show the M^g distribution at the 
benchmark point. The peak of this distribution shows 
a well defined peak for 50 fb^^ luminosity. 




Meff (GeV) 

FIG. 5: Distribution of M^s at a stop coannihilation bench- 
mark point. The distribution is fitted with a Gaussian func- 
tion to find the peak. The value of the AI^^^^ observable is 
1073.01 ± 8.72(Stat.) GeV. The luminosity is 50 fb"^ 

In Table [V] we show the end points of Mtt and Mjtt 
distributions, peak of the Afcff distribution, and the 
slopes of the pr distributions for the taus for our bench- 
mark point. The statistical uncertainties range between 
0.4% - 3.8%. 

TABLE V: Kinematic observables at 50 fb^^ for a stop coan- 
nihilation benchmark point. All masses are in GeV. 



Observable 


Value 


50 fb~' Stat. 


100 fb"' Stat. 


slope{pTMsh) 
slope{pT.iow) 

en 


49.71 
269.09 
-0.0522 
-0.1178 
1073.01 


±0.20 
±3.18 
±0.0021 
±0.0040 
±8.72 


±0.14 
±2.25 
±0.0014 
±0.0028 
±6.17 



C. Observables for the Determination of Third 
Generation Squark Masses: M§w and Mj^ 

To probe the third generation squark masses we need 
to involve b quarks. The relevant decay chain associated 
with the dominant production process for the reconstruc- 
tion of third generation squarks in the stop coannihilation 
region is 

g — >b + b — >i+W + b — > Xi + c + b + W. (13) 

These signals are characterized by high px jets accom- 
panied by a and pr- 



We will be constructing the distributions M^w and 

The cuts for the analysis are 

{i) pT > 180 GeV; 

(ii) Number of jets: iVjct > 4; 

{in) Leading jet cuts: The first two leading jets each 
have pt > 200 GeV in \?]\ < 2.5. They could be gluon, 
light-flavor, or b jets; 

{iv) Soft jet cuts: Any jets with visible px > 30 GeV 
in 1 77 1 < 2.5 are accepted in the analysis. This includes 
&-tagged jets; 

(f ) PTjetl + PT,jet2 + pT ^ GOO GcV; 

(vi) For Mbw, at least one tight b jet is required. 

The flrst step in the analysis is the reconstruction of 
the W boson. The W appears in the detector as two jets 
whose invariant mass falls in the W mass window (65 
GeV < Mjj < 90 GeV). We thus choose soft jet pairs 
(from the third leading jet and below) which are not b- 
tagged, with 0.4 < AR < 1.5. The jets are put into two 
categories: those which are manifestly in the W window, 
and those that fall within the sideband window (40 GeV 
< Mjj < 55 GeV or 100 GeV < Mjj < 115 GeV). BEST 
is then performed for the two categories, to get rid of 
uncorrelated jet background. After this, the sideband 
subtraction is performed to obtain the W mass. 

Once the W is reconstructed, it is paired up with jets 
to form the Mi,w and Mjw distributions. 
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FIG. 6: Distribution of Mjw at a stop coannihilation bench- 
mark point. W is first reconstructed with two jets whose 
invariant mass falls in the W window. The reconstructed W 
is then combined with a non 6-tagged soft jet of rank three 
or lower from the same event, to produce the same-event blue 
(grey) histogram. The W is combined with a soft jet from a 
difi'erent event to produce the bi-event filled dot-dashed blue 
(grey) histogram, which is normalised to the shape of the 
long tail of the same-event histogram. The same-event minus 
bi-event subtraction (BEST) produces the black subtracted 
histogram. The subtracted histogram is fitted with a straight 
line to obtain the endpoint. The result for the endpoint is 
287.55 ± 0.74(Stat.) GeV. The luminosity is 50 fb"^ 

For the Mjw distribution, we pair the W with a non 
6-tagged soft jet, whose rank is three or lower. This is be- 
cause we are in the stop coannihilation region. In Figure 
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FIG. 7: Distribution of Mtw at a stop coannihilation bench- 
mark point. W is first reconstructed with two jets whose in- 
variant mass falls in the W window. The reconstructed W is 
then combined with a 6 jet of any rank from the current event, 
to produce the same-event pink (grey) histogram. Events with 
Mbw < 200 GeV are discarded to remove the top peak. The 
W is combined with a b jet from a different event to pro- 
duce the bi-event filled dot-dashed pink histogram, which is 
normalised to the shape of the long tail of the same-event his- 
togram. The same-event minus bi-event subtraction (BEST) 
produces the black subtracted histogram. The subtracted his- 
togram is fitted with a straight line to obtain the endpoint. 
The resuh for the endpoint is 325.67 ± 4. 50(Stat.) GeV. The 
luminosity is 50 fb~^. 



l6j we show the Mjw distribution at the benchmark point, 
finding a well defined end-point for 50 fb~^ luminosity. 

For the Mhw distribution, we pair the W with a b jet 
of any rank from the current event. Note that the rele- 
vant 6 jet required to construct this observable need not 
be a leading jet; in fact, the leading jet will typically be 
non ^-tagged. After pairing the W with the b jet, we do a 
further BEST to get rid of uncorrelated b jets. This gives 
the final signal for M^vi/- However, the b+W signal shows 
the presence of the unwanted top peak, which comes from 
t — > b + W . The top window is removed from the fi- 
nal signal, by discarding events with M^w < 200 GeV. 
In Figure [7] we show the Mh^r distribution obtained at 
the benchmark point, finding the end-point for 50 fb~^ 
luminosity. 

In Table |VH we show the endpoint values obtained 
from these distributions. The statistical uncertainties 
range between 0.2% — 1.4%. The statistical uncertainty 
is larger for Mtw due to the b jet. 



TABLE VI: Kinematical observables M'. 



end 
jW 



fb ^ for a stop coannihilation benchmark point. All masses 
are in GeV. 



Observable 


Value 


50 fb"' Stat. 


100 fb"' Stat. 


IVljW 


287.55 
325.67 


0.74 
4.50 


0.52 
3.18 
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IV. DETERMINATION OF PARTICLE MASSES, 
GAUGINO UNIFICATION SCALE, AND UV 
MODEL PARAMETERS IN THE STOP 
COANNIHILATION REGION 



The observables measured in Section IIIII are used to 
determine the masses of the gluino, Xi: X2' ^^"^ 9' 
which lead to the determination of model parameters 
(a, 7713/2, tan/3, rim, fi_f/)- We will also determine the b 
and i masses later in this section. 

Theoretically, the functional dependences of the ob- 
servables in Section [lTT| are expected to be: 

slope{pT, am) = PTAM{m^o,m^o) 
slope{pT,dia) = PT,dis{mf,m^o,m^o) 
{mq^,mg) 



11 rpcak 11 rpcak 



(14) 



The masses are varied independently around the 
benchmark point, and the full simulation of the collider 
experiment and determination of the observables is per- 
formed each time. Thus, the functional dependence of 
the observables on the masses is determined, along with 
their uncertainties. 

This set of equations can be used to solve for the 
masses. In general, the dependence of the observables on 
the masses is expected to be non-linear, corresponding 
to multiple solutions of the masses. However, near the 
benchmark point, we found that for most masses, linear 
functions fitted the dependence quite well. We obtained 
two solutions only for the stau mass, and one was quite 
far from the benchmark point and thus rejected. We will 
have more to sa y on these issues in Section |V A[ 

we show the solution to the g,qL,T,X2, 
^ luminosity. We also show the 



VII 



In Table 
and Xi masses for 50 fb 
statistical uncertainties, which range between 1% — 6.6%. 



TABLE VII: Solution to masses at 50 fb ^ for a stop coanni- 
hilation benchmark point. All masses are in GeV. 



the full model. As we have stressed in Section|Tj the point 
of this work is to show that the contours of the media- 
tion scheme (here, mirage pattern) show up unmistak- 
ably from the determination of certain masses. We will 



demonstrate this in Section IV C where we will show the 



gaugino unification scale. However, a set of masses may 
point back at different models, and in particular, we will 
not attempt to perform a distinguishability study of our 
particular model. When we solve the model parameters 
for our particular example (KKLT), our methods will re- 
flect this - we will start from the model parameters, flt 
with the spectrum, and zero in on the masses we solved. 

We use the Nelder Mead method, a commonly 
used nonlinear optimization technique, to do so. An 
initial simplex is chosen on the parameter space 
{777.3/2, a, rim, nn, tan/3, } and ISAJET [17] is used to flnd 
the corresponding mass. The scan over parameter space 
proceeds until the masses calculated by ISAJET are close 
to the solved values. 



Table VIII shows the corresponding model parameters 
we determined from 777^, mqj^, m^o, rrif-^ and 777^^0 for 50 
fb~^ and 100 fb~^ luminosities. At 50 fb""'^, the statistical 
uncertainties range between 5.6% — 17%. 



TABLE VIII: Solution to model parameters at a stop coan- 
nihilation benchmark point at 50 fb~^. Masses are in GeV. 



Parameter Value 


50 fb"" Stat. 


100 fb"^ Stat. 


a 


4.58 


±0.21 


±0.14 


m-i/2 


13717 


±688 


±517 




0.106 


±0.015 


±0.015 


riH 


0.578 


±0.095 


±0.091 


tan/3 


28.76 


±1.65 


±1.36 



Using the model parameters from Table |VIIH we cal- 
culate the relic density and we find: 

= 0.096 ± 0.029 . (15) 

The relic density is determined to an accuracy of 31%. 



Particle Mass 50 fb"^ Stat. 


100 fb"' Stat. 


g 646 


-14,±19 


-11,+14 


QL 638 


-34,±42 


-23,+39 


f 318 


-3,±3 


-3,+3 


X2 333 


-7,+ll 


-6,+8 


X? 276 


-8,+13 


-7,±10 



A. Determination of UV Model Parameters and 
Relic Density 

Having determined the masses of g, x?, X21 ^^'^ 9' 
can use these to numerically solve for the parameters of 



B. Determination of b and t Masses 



The observables and M^^ are used to deter- 

mine the third generation squark masses, once the other 
masses have been obtained with the observables previ- 
ously determined. 

Theoretically, the functional dependences are Mi,w = 
Mi,iY{mf^,m^,mg) and Mj^r = Mjw{mi,m^,m^o). As 
before, the masses are varied independently around the 
benchmark point, and the collider experiment and deter- 
mination of M^^ and is performed each time. 

We show the masses of the third generation squarks 
with uncertainties in Table HXl The statistical uncertain- 
ties range between 2.5% — 11.3%. 
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TABLE IX: Solution to stop and sbottom masses at 50 fb ^ 
for a stop coannihilation benchmark point. All masses are in 
GeV. 



Particle Mass 50 fb"^ Stat. 


100 fb" 


^ Stat. 


b 531 -60, +60 


-47, 


+47 


i 326 -5, +8 


-4, 


+7 



C. Gaugino Unification 

Upto this point, we have displayed our techniques of re- 
constructing masses at the benchmark point given in Ta- 
ble [TTj Our techniques work perfectly well at benchmark 
points with higher gluino mass, as preferred by current 
LHC data. Higher luminosity is of course required to 
obtain endpoints. Below, we choose such a benchmark 
point with nig ^ 1.2 TeV, reconstruct gaugino masses, 
and show the gaugino unification scale. 



TABLE X: Model parameters and spectrum at a new stop 
coannihilation benchmark point with heavier gluino. All 
masses are in GeV. 



Parameter Value Particle Mass 



a 


3.8 


9 1187 


mz/2 


34800 


X2 740 


Tim 


0.0 


X? 666 


riH 


0.5 


f 721 


tan/? 


28 


q 1189 



As mentioned in Section [Tj gaugino masses may be ob- 
tained at lower luminosity, as we elaborate below. 

The solution of gaugino masses typically requires the 
values of four observables (pT,diff ,Pt,am, and 
M™'*), as is clear from the solutions we displayed in Table 



VII One could have used the endpoint of the invariant 
mass distribution Mjr of the leading jet and the more 
energetic r as an observable. Our analysis of this ob- 
servable indicates that one requires more luminosity to 
find its endpoint, relative to the luminosity required to 
determine the observables pt, am and pT,diff, because of 
the associated jet subtractions. 

Since for the benchmark point the r's are quite close 
in mass, we use pt.am and PT,diff, obviating the need 
to use observables which involve jet information. Then, 
XiJ X2j and f masses can be obtained from PT.AM,PT.diff , 
and M™*^ observables. Since these observables do not 
involve jet subtractions, they give acceptable endpoints 
at lower luminosity. On the other hand, M^^'^^ roughly 
gives the gluino mass. 

In Table pClH we show the values of the gaugino masses 
at 50 fb~^ for the new benchmark point. The statistical 
uncertainties range between 1.2% — 4.2%. 

In Figure [Sj we show the running of the gauginos and 
the gaugino unification for our benchmark point. Clearly, 
gaugino unification occurs below the GUT scale, estab- 
lishing the footprints of anomaly contributions to super- 
symmetry breaking. 



TABLE XII: Solution to gaugino masses at 50 fb~^ for the 
new stop coannihilation benchmark point with heavier gluino. 
All masses are in GeV. 



For this benchmark point, a luminosity of 200 fb"^ is Particle Mass 50fb ^ Stat. 

required to solve for all the masses, following the tech- 
niques we have shown in this paper. We show the masses 



we obtained for this benchmark point in Table XI The 



statistical uncertainties range between 0.9% — 14.7%. 



9 


1181 


± 50 




738 


± 15 


X? 


649 


± 20 



TABLE XI: Solution to masses at 200 fb~^ for a stop coanni- 
hilation benchmark point with heavier gluino. All masses are 
in GeV. 



Particle Mass 


Stat. 


t 690 


± 6 


b 1002 


± 126 


T 717 


± 10 


q 1133 


-132, +167 



The model parameters are solved as before, using the 
masses above as well as the gaugino masses solved at 200 
fb^^. The solution to the model parameters at 200 fb^^ 
for the new benchmark point are a = 3.84 + 0.11, 7713/2 — 
34463±805 GeV,7i™ = 0.078±0.007, 77^ = 0.592+0.036, 
and tan/3 = 27.2 + 1.94. The relic density is found to be 
VtK^ =0.23 + 0.13. 



The theoretically expected gaugino unification scale, 
for this benchmark point, is given by 

/777 \2-47/a 

Munif = Mgut y-^j - 10^ GeV . (16) 

Note that the expression differs from the ones in [8] by 
the rescaling in Eq. [8) 

D. Systematic Errors 

We also study an impact due to a ±3% uncertainty 
on the energy scale for jets, taus, and missing transverse 
energy to estimate the systematic uncertainties indepen- 
dently of luminosity. 

Table |XIII| shows the percentage systematic errors in 
the determination of the observables and masses. We 
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FIG. 8: Running of gaugino masses. The vertical axis is 
masses in GeV. The luminosity is 50 fb~^. The shaded lines 
from the top down represent, at the left edge of the graph, 
the gluino, wino and bino masses respectively. 



find that the systematic errors for the observables are be- 
tween 1% — 6%, which is larger than the statistical error 
0.2% — 3% at 50 fb""'^. The systematic errors of masses, 
which vary from 2% — 19%, are also large compared to the 
statistical errors, which range between 1% — 11% at 50 
fb~^. Since our assumption of ±3% on the energy scale 
could be changed in actual experimental conditions, we 



just show the size of errors in Table XIII as a reference 
and don't use them in the analyses performed in the pre- 
vious sub-sections. 



point with heavier gluino mass preferred by current LHC 
data. 

The distributions used in the stau coannihilation re- 
gion are Mrr, Mjrr, Mjr, pr.low, and M^s- Note that 
we only use the pr.iow corresponding to the lower energy 
T from f — > Xi + since the higher energy r doesn't 
show enough variation as the masses are varied. 

We use Mjr , which is formed by combining the higher 
energy r with a leading jet from the same event. OS— LS 
and BEST are performed as usual. Theoretically, one 
expects M:^f = Mjf {niq^ , m^o , ruf) . 

In Figures |9] [TOj |lT] [12] and we present the distri- 
butions obtained. 



180 
> 160 




40 60 80 100 120 140 

Mx. (GeV) 



FIG. 9: Distribution of Mtt at a stau coannihilation bench- 
mark point. The notation is the same as in Figure [l] The 
endpoint obtained is 90.70 ±0.54(Stat.) GeV. The luminosity 
is 100 fb-^ 



TABLE XIII: Systematic percentage errors in observables and 
masses for a benchmark point in the stop coannihilation re- 
gion. 



Observable Error (%) Mass Error (%) Parameter Error (%) 





± 3.0 


m-g 


± 3.9 


Q 


± 5.9 




± 4.5 




-2.4,-^4.2 




± 5.9 


s;ope(pT,high) 


± 5.2 


m^o 


± 4 




± 12.5 


slope{pT,iov,) 


± 1.0 


rUf 


± 2.6 


riH 


± 17 


^peak 
eiT 


± 1.3 


ruq 


-13,4-15 


tan/3 


± 6.1 




± 5.9 




-2.5,-1-3.7 








± 3.7 




± 19 
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V. ANALYSIS IN THE STAU 
COANNIHILATION REGION 

In this Section, we determine masses and model pa- 
rameters for a different point in parameter space, where 
only stau coannihilation is operational. We first show our 
analysis for the benchmark model point shown in Table 
|III| In Section |V C[ we will also analyse a benchmark 



FIG. 10: Distribution of Mjrr at a stau coannihilation bench- 
mark point. The notation is the same as in Figure |2] The 
endpoint obtained is 479.53±3.45(Stat.) GeV. The luminosity 
is 100 fb-^ 

In Table |XIV( we present the endpoints obtained from 
the Mtt, MjTT, and Mjr distributions, the slope of the 
visible pt distribution of the lower energy r, and the peak 
of the Mcff distribution, at 100 fb"'^ for the stau coanni- 
hilation benchmark point. The statistical uncertainties 
range between 0.6% - 4.7%. 
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FIG. 11: Distribution of pT,iow at a stau coannihilation bench- 
mark point. The notation is the same as in Figure |4] The 
slope obtained is -0.0849 ± 0.0041(Stat.). The luminosity is 
100 fb^^ 
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FIG. 12: Distribution of M^fi at a stau coannihilation bench- 
mark point. The notation is the same as in Figure |5] The 
peak obtained is 1257.26 ± 10.33(Stat.) GeV. The luminosity 
is 100 fb-\ 



TABLE XIV: Kinematical observables at 100 fb"^ for the 
Stau coannihilation benchmark point. All masses in GeV. 

Observable Value 100 fb"^ Stat. 



slope{pT,r) 



90.70 
479.53 
-0.0849 
1257.26 
448.40 



±0.54 
±3.45 
±0.0041 
±10.33 
±16.20 



A. Masses and Gaugino Unification 

After obtaining the values of the observables at the 
benchmark point as quoted above, we varied the masses 
around the benchmark point. For each mass selection, 
we simulated the experiment and obtained all the ob- 
servables. This gave us the observables as a function of 
the masses. We then iteratively solved for the masses. 



FIG. 13: Distribution of iVf,> in jrr system. The blue (grey) 
histogram is obtained by combining the higher energy r of 
the OS— LS r pair with a leading jet of the same event. The 
filled dot-dashed histogram is the bi-event histogram, con- 
structed by combining the aforesaid t with a leading jet from 
a different event and is normalised to the long tail of the same- 
event histogram. The same-event minus bi-event subtraction 
(BEST) produces the black subtracted histogram. The sub- 
tracted histogram is fitted with a straight line to obtain the 
endpoint. The resuh for the endpoint is 448.40 ± 16. 20(Stat.) 
GeV. The luminosity is 100 fb"^ 



Obviously, one should expect more than one set of mass 
solutions in this method, since the dependence of the 
observables on the masses near the benchmark point is 
non-linear. In this case, we found two sets of solutions 
with quite disparate values of masses. 

This is a general problem, and we mention some ar- 
guments to reject the "wrong" solutions. Firstly, one 
set of solutions had masses far away from the bench- 
mark point, where the expansion was not performed to 
begin with. Secondly, this set had a much larger stau- 
neutralino mass difference, which did not satisfy the relic 
density constraint. 

Finally, one can reject masses by appealing to the 
UV model and determining more observables, although 
this goes somewhat against the bottom-up philosophy we 
have pursued in this paper. Thus, we can use this set of 
mass solutions and solve for the UV model parameters, 
which we can then use to solve the spectrum of the model, 
in particular obtaining the b quark mass. Given the spec- 
trum, we can obtain the theoretical value of a particular 
observable, say M^Wj and on the other hand, measure 
it using the techniques outlined in the paper. For the 
incorrect mass solution set, the values thus obtained will 
be very different. 

In any case, after rejecting the first set of solutions, the 
second set itself had two sets of values of x? and f masses, 
which were clos e to each other. One set of solutions is 
XV The other solutions to the Xi and 
fo = 294 GeV and mf 



shown in Table 
f masses are m 



337 GeV. All 

combinations of masses (four in all) are used to solve the 
model parameters, as we outline in the next subsection. 
We display the solution set in Table |XV[ obtained at 
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luminosity 100 fb ^. The statistical uncertainties are 
2.7% - 6.4%. 

TABLE XV: Solution to masses at 100 fb~^ for the stau coan- 
nihilation benchmark point. All masses are in GeV. 



Particle Mass 100 fb"^ Stat. 

I 895 -35, +50 

qL 845 -36, +24 

X2 388 -9, +25 

f 298 -8, +8 

X? 274 -10, +10 



C. Gaugino Unification 

Although upto this point, we have displayed our results 
in the stau coannihilation region at the benchmark point 
given in Table |IV[ our methods are applicable at other 
benchmark points with higher gluino mass, preferred by 
current LHC data. We display such a benchmark point 
for stau coannihilation below, and proceed to show gaug- 
ino unification. 

TABLE XVII: Model parameters and spectrum at a new 
stau coannihilation benchmark point with heavier gluino. All 
masses are in GeV. 



Parameter Value Particle Mass 



B. Determination of UV Model Parameters and 
Relic Density 

Having determined the masses of g, XitX2i "^i q, we 
can use these to numerically solve for the parameters of 
the full model. 

We use the Nelder-Mead method to do so. An 
initial simplex is chosen on the parameter space 
{™3/2j Q^j n-m, tan/3, } and ISAJET is used to find the 
corresponding mass. The scan over parameter space pro- 
ceeds until the masses calculated by ISAJET are close to 
the solved values. 

In the previous subsection, we have discussed two solu- 
tions each for Xi and f . We solved the model parameters 
for each combination of masses (four in all) , and took the 
average of the values of the model parameters thus ob- 



tained. These average values are displayed in Table pCVT 
The statistical uncertainties are 9% — 40%. 



TABLE XVI: Solution to model parameters at the stau coan- 
nihilation benchmark point. For each mass solution, the 
model parameters are solved, and the average of all the dif- 
ferent sets of solutions is shown in the table. Masses are in 
GeV. 



Parameter Value 


Stat. 


a 


7.42 


± 0.58 




10171 


± 882 




0.52 


± 0.09 


riH 


1.17 


-0.07, +0.22 


tan/? 


33.1 


± 7.8 



Using the model parameters from Table |XVlj we cal- 
culate the relic density and we find: 



= 0.17t°j^ 



(17) 



The relic density is thus determined to an accuracy of 
59%. 



a 


10 


9 


1183 




9500 


xl 


499 


Tim 


0.5 


X? 


337 


riH 


1.0 


f 


361 


tan/3 


37 


q 


1119 



At this benchmark point, we found acceptable end- 
points at 250 fb~^, and solved for the observables and 
masses. The masses are mf = 346 ± 12 GeV and 
rriq — 1138 ± 72 GeV. The masses of the gauginos 
can also be determined, although we show their val- 
ues at a lower luminosity below. After determining 
the masses, the observables can be solved and they 
are a = 9.95 ± 0.99, = 9636 ± 910 GeV, = 

0.56 ±0.09, UH = l-Wlo.os, and tan/3 = 40.0+7.0. The 
relic density for the heavy gluino point is found to be 

nK^ = omtall- 

The gaugino masses can be obtained at lower luminos- 
ity, as spelled out in Section |lj 

In the stau coannihilation region, one can use the prox- 
imity of the stau and lightest neutralino masses to solve 
gaugino masses using fewer observables, in particular 
avoiding the observable M™'*, which typically requires 
higher luminosity due to jet subtractions. This is remi- 
niscent of what we did in the previous sections with stop 
coannihilation, although there we used the proximity of 
r pt's. 

We take approximately 



TTif ^ m^o + 15 



Thus, assuming that A'l^^'^^ is approximately determined 
by m^, it is clear from Eq. 14 that one can solve for 
mg,m^o,m^o from M^^^^ ,pT,iow and M^t. 

These observables have acceptable endpoints and 
peaks down to 15 fb~^. The main difference from the 
full set of observables is the Mj^ distribution, which re- 
quires a higher luminosity to obtain acceptable endpoints 
due to jet subtractions. 

In Table [XVinj we show the gaugino masses obtained 
at 15 fb~^ for the new stau coannihilation benchmark 
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TABLE XVIII: Solution to gaugino masses at 15 fb"^ for the 
new stau coannihilation benchmark point with heavier gluino. 
AH masses are in GeV. 



Particle Mass 15fb"^ Stat 



g 

X? 



1186 
527 
317 



± 84 
-30, +60 
± 33 



point with heavier gluino. 
are 7% - 11.4%. 



The statistical uncertainties 



In Figure 14 we show gaugino unification for the new 
benchmark point at 15 fb^^. 
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FIG. 14: Running of gaugino masses. The vertical axis is 
masses in GeV. The luminosity is 15 fb~^. The shaded lines 
from the top down represent, at the left edge of the graph, 
the gluino, wino and bino masses respectively. 



VI. CONCLUSIONS 

If supersymmetry is discovered at the LHC, the next 
step would be to understand as much as possible about 
the underlying model of supersymmetry breaking and 
mediation. How should one proceed? 

One option is to start from specific models defined by 
a tractable number of parameters at high energy, and 
fit to observables measured in the collider. However, we 
have stressed that a given set of measured observables 
will generally point back at different UV models, dis- 
tinguishing among which may be challenging and even 
intractable (depending on the complexity of the model 
and the number of reliable observables measured). 

We have stressed on a bottom-up approach to this 
question: first identify sectors in the low energy spec- 
trum which carry the imprints of mediation schemes in 
a relatively model-independent manner, next construct 



observables which carry information about these masses, 
and finally using the observables solve for the masses and 
identify a mediation scheme. 

The gaugino sector offers a particularly fruitful sector 
in this program, since the determination of this sector 
can unambiguously show the imprints of anomaly contri- 
butions to supersymmetry breaking. We therefore take 
an example of a model where anomaly contributions are 
significant (supersymmetry breaking in KKLT compact- 
ification) and proceed to obtain the gaugino masses, by 
constructing suitable observables. 

We rely on the construction of kinematic observables 
arising from cascade decays using endpoint techniques, 
since these are a particularly sharp tool for such diagno- 
sis. 

While the gauginos are important in the paradigm out- 
lined above, they form only a subset of the masses we de- 
termined in this paper. Given a neutralino dark matter, 
the determination of other masses becomes a necessity 
when one wants to establish mechanisms to satisfy the 
relic density. This is perfectly amenable to the kind of 
bottom-up approach we have taken: in our examples, we 
needed to establish that the stop-neutralino and stau- 
neutralino coannihilation mechanisms were operational. 
To do so demanded the determination of stop and stau 
masses, in addition to the neutralino. 

We summarize our work below: 

(z) We constructed observables using different combi- 
nations of the jets, r's and Ws in the final states. We 
constructed these observables in the stop-neutralino and 
stau-neutralino coannihilation regions. 

In the stop coannihilation region, we constructed two 
new observables, MJJ^ and M^^, to deter mine the 

masses of the i and the h. The determination of the i 
and h masses is especially challenging, since both decay 
to h quarks. To make the problem worse, stop decays 
to missing energy by emitting a low energy jet in this 
coannihilation region. We showed how to reconstruct the 
stop mass in order to establish the fact that we are in the 
stop-neutralino coannihilation region. 

We also constructed two other new observables, pt,am 
and pT.diff which are combinations of the pr of two op- 
posite sign r's. These observables are particularly useful 
when the sample shows the presence of two opposite sign 
r's that are close in mass. 

[ii) Using the observables, we solved for the masses of 
the gluino, the two lighter neutralinos, squarks (of the 
first two generations) and the lightest stau. We also de- 
termined the lighter stop and sbottom masses in the stop 
coannihilation region. 

{Hi) We used the neutralino and gluino masses to de- 
termine the gaugino unification scale and discern the ef- 
fects of anomaly mediation from a bottom-up approach, 
as elaborated above. 

{iv) Using all the masses, we determined (a) the model 
parameters and (b) tested whether we are in a coanihi- 
lation region and dark matter relic density is satisfied. 
We find that relic density can be determined with an 
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accuracy of 31% in the stop and 59% in the stau coanni- 
hilation region. 

As far as the diagnosis of the supersymmetry break- 
ing mediation scheme is concerned, the most important 
future direction is obviously to expand our techniques 
to parts of the parameter space away from the coanni- 
hilation regions. Our observables enabled us to solve for 
the lightest neutralinos, and for the statements on medi- 
ation scheme that we made above, they were assumed to 
be gauginos (this is guaranteed in the coannhilation re- 
gions) . It would be very interesting to test the paradigm 
outlined in this paper at regions where the Higgsino com- 
ponent is not insignificant. Also interesting would be to 
find bottom-up techniques to further distinguish between 
different schemes that correspond to the same gaugino 
mass pattern. We leave these questions for future work. 
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